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LONG-TERM  GOALS 

The  overall  long-term  goal  for  this  project  is  to  determine  if  and  how  acoustic  coherent  backscatter 
enhancement  (ACBE)  can  be  used  for  classification  of  active  sonar  returns  in  a  wide  variety  of  ocean 
environments.  During  its  first  two  years,  this  project  has  focused  on  simulations  of  acoustic  multiple 
scattering  from  aggregations  of  omni-directional  point  scatterers  to  determine  the  parametric  realms  in 
which  ACBE  might  be  observed,  and  its  characteristics  when  it  is  observed. 

OBJECTIVES 

The  detailed  objectives  of  the  current  research  effort  are  to  determine  the  parametric  dependence  of 
ACBE  peak  amplitude,  peak  emergence  rate  as  the  number  of  observations  increases,  peak  angular 
width,  and  peak  time  dependence.  Here  the  independent  parameters  are  the  range  between  the 
scattering  aggregation  and  the  receiving  array,  the  receiving  array  characteristics,  incident  wave 
characteristics  (wave  type,  waveform,  frequency,  bandwidth,  duration),  and  aggregation  characteristics 
(scatterer  cross  section  and  mean  spacing,  overall  aggregation  size  and  shape,  etc.).  Eventually, 
underwater  waveguide  characteristics  will  be  considered  as  well. 

APPROACH 

The  current  approach  involves  numerical  evaluation  of  the  fundamental  equations  of  multiple 
scattering  from  an  aggregation  of  omni-directional  point  scatterers1.  If  yAj)  is  the  harmonic  acoustic 
pressure  field  at  frequency  co  at  the  point  r  and  y/o(r)  is  the  harmonic  field  incident  on  the  aggregation 
of  scatterers  located  at  r„,  then 

N 

V<r)=  ^0(r)+  ^(r)  =  ^0(r)+Z«IVf.(r.)G(r»r«)  >  C1) 

n  1 

where  y/[jr)  is  the  scattered  field  and  is  given  by  the  sum  in  (1),  A  is  the  number  of  scatterers,  g„  is  the 
scattering  coefficient  of  the  nth  scatterer,  i//„(r„)  is  the  field  incident  on  the  nth  scatterer, 

V'(r„)  =  n(r„)+  Z  (2) 

J 

and  G(rn,  rj)  is  the  free-space  Green's  function  between  the  locations  r„  and  r y, 
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where  ko  =  cole  is  the  wave  number  magnitude  of  the  incident  field,  c  is  the  sound  speed,  and  i  -  a/-1. 
When  the  incident  field  and  the  scattering  coefficients  are  known,  (2)  can  be  written  N  times,  once  for 
each  scatterer  1  <n<N,  and  these  N  algebraic  equations  can  be  solved  to  determine  y/n( r„).  The  total 
field  at  any  location  is  then  recovered  from  (1)  and  (3)  using  the  known  y/o(r),  the  known  g„,  and  the 
calculated  y/n(r „).  This  formulation  is  akin  to  the  direct  boundary-integral  formulation  of  computational 
acoustics  with  one  computational  element  assigned  to  each  scatterer.  The  computational  burden  of  this 
approach  is  set  by  the  inversion  of  the  fully-populated  N-by-N  algebraic  system  that  determines  y/„(rn). 
In  the  current  investigation,  the  scatterers  are  placed  randomly  with  an  average  spacing  s,  and  are 
identical  without  internal  lossless  so  the  N  scattering  coefficients  are  all  the  same 
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where  crs  is  the  scattering  cross  section.  In  the  current  investigations,  crv  is  considered  an  independent 
parameter  within  the  constraint  imposed  by  conservation  of  acoustic  energy:  <r,  <\ni k%  =  X1  n.  For 
the  current  calculations,  the  incident  field  y/o(v)  is  a  plane  wave  with  amplitude  A  and  wave  number 


vector  k0,  and  the  backscatter  direction  (</>=  <p  =  0)  is  defined  with  respect  to  kQ  as  shown  in  Figure  1. 

Rectangular  aggregations  of  scatterers  with  width  and  depth  dimensions  h  and  w  have  been  the  primary 
geometry  considered  this  fiscal  year  because  this  geometry  readily  allows  both  near-field  and  far-field 
investigations  with  a  linear  receiving  array  of  length  L  when  the  number  of  scatterers  is  limited  to  a 
few  thousand. 


To  search  for  the  presence  or  absence  of  ACBE,  the  scattered  field  ^(r)  predicted  by  (1)  is  calculated 
at  the  elements  of  the  receiving  array,  and  aggregation-array  distance  R  is  varied  to  put  the  receiving 
array  in  the  near-field,  where  the  array's  beam-steering  angle  (f>  is  relevant,  or  in  the  far-field  of  the 
aggregation,  where  the  azimuthal  scattering  angle  <p  is  relevant.  For  near-field  calculations,  the  results 
are  presented  as  B((/))l[B\lve  vs.  (f>,  where  B(<f>)  is  the  beamformed  output  of  the  receiving  array,  and 
[B]ave  is  the  array's  average  beamformed  output  in  directions  near  backscatter  but  excluding  the  ACBE 

peak.  For  far-field  calculations,  the  results  are  presented  as  \yy(jp)\2  | yyj2  vs.  <p,  where  \\f/ |2  is  the 

average  mean-square  pressure  at  that  aggregation-receiver  range  in  directions  near  backscatter  but 
excluding  the  ACBE  peak.  The  independent  parameters  of  these  investigations  are  A,  ko,  crs,  s,  R,  L,  <f> 
or  <p,  w,  and  h.  Interestingly,  the  current  problem  can  be  stated  in  dimensionless  terms  by  normalizing 
field  values  with  the  incident  wave  amplitude  A,  and  by  using  wave  number-scaled  lengths.  This 
reduces  the  number  of  independent  parameters  on  which  ACBE  may  depend  from  nine  to  seven. 

The  current  focus  of  this  research  effort  is  on  determining  and  understanding  the  parametric 
dependence  of  the  ACBE  peak's  height  and  width  on  the  seven  dimensionless  parameters  (kQ  Sjo\.,  kos, 

pR,  pL,  (f>  or  <p,  kow,  and  k(,h),  a  task  that  has  proved  challenging.  In  particular,  the  current  simulations 
show  the  ACBE  peak  enhancement  may  be  10  dB  (or  more)  above  the  scattered  intensity  in  other 

directions  depending  on  the  values  of  koyjo\.  and  ps.  A  maximum  enhancement  of  3  dB  was  expected 

based  on  prior  theory2.  Thus,  understanding  ACBE  peak  amplitude  dependence  on  kQ^Jas,  kos  and 

aggregation  shape  has  been  the  focus  of  this  year's  work.  Investigations  to  date  involving  the  remaining 
geometric  parameters  suggest  they  play  more  easily  understood  roles. 
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Figure  1.  Geometrical  configuration  of  the  multiple  scattering  simulations.  A  plane  wave  with 
amplitude  A  and  wave  number  vector  k0  impinges  on  a  rectangular  aggregation  of  scatterers  with 

dimensions  w  and  h.  Here  (pis  the  scattering  angle  defined  from  the  origin  of  coordinates,  and  <j)is 
the  beam-steering  angle  defined  from  the  broadside  direction  of  the  receiving  array  of  length  L  that 
lies  a  distance  R  from  the  aggregation.  The  backscatter  direction  is  (p=  <f>=  0. 


These  ACBE  investigations  are  the  current  doctoral  research  of  Ms.  Adaleena  Mookerjee.  She  is  a  US 
Citizen  and  a  Ph.D.  candidate. 

In  addition,  an  undergraduate  team  supported  by  NAVSEA  through  the  Naval  Engineering  Education 
Center  (NEEC)  has  developed  a  tank  experiment  (42"  diameter,  42"  depth,  16  receiving  hydrophones) 
that  may  eventually  be  suitable  for  short-range  ACBE  experiments  in  the  frequency  range  from  30  kHz 
to  120  kHz. 

WORK  COMPLETED 

The  simulation  code  was  subjected  to  three  validation  tests  involving  three  different  geometries.  First 
acoustic  energy  conservation  was  checked  for  scattering  from  an  aggregation  of  1020  scatterers  placed 
randomly  within  a  sphere  approximately  five-wavelengths  in  diameter.  For  these  simulations,  the 

scatterers  were  strong,  k0Sjo\.  =  3.5,  and  relatively  close  to  each  other,  k(,s  =  2.5.  When  the  incident 

plane  wave  is  included  the  net  acoustic  power  crossing  a  surface  that  encloses  any  aggregation  of 
scatterers  should  be  zero.  The  net  power  recovered  from  the  simulations  was  0.6%  of  the  power 
scattered  by  one  scatterer.  Thus,  the  simulations  are  very  nearly  energy  conserving. 
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Next,  a  linear  arrangement  of  13  evenly  spaced  scatterers  was  illuminated  from  the  end-fire  direction 
and  the  angles  of  maximum  scattered  intensity  in  the  far-field  of  this  linear  aggregation  were 
determined  when  kos  =  8.38  and  16.8.  The  angles  of  increased  scattered  intensity  recovered  from  the 
simulations  successfully  matched  those  determined  from  Bragg  scattering  formulae.  Thus,  the  relative 
phasing  of  the  scattered  waves  in  the  simulations  appears  to  be  correct. 

The  remaining  simulation  test  involved  comparing  the  simulation  results  with  theory  and  experiments, 
and  is  described  in  the  next  section. 

RESULTS 

Sample  far-field  ACBE  simulation  results  are  shown  in  Figure  2  for  a  cubical  aggregation  with  a  five- 
wavelength  edge,  256  to  1024  Monte -Carlo  trials,  k0Sjo\.  =  3.5,  and  four  different  average  scatterer 

/ 1  ,2  \  f  /\  |2  \ 

spacings:  kos  =  11.7,  5.5,  3.2,  and  2.5.  Here,  the  vertical  axis  is  y/s(R,<p)  yfjjt, <p)  and  the 

\  /  /  \  /  ave 

horizontal  axis  is  the  far-field  scattering  angle  (p.  The  (  ^-brackets  denote  an  ensemble  average  over 

Monte-Carlo  trials.  The  central  peak  at  cp  =  0  for  k(,s  =  3.2  and  2.5  is  produced  by  coherent 
backscattering  enhancement.  Interestingly,  the  enhancements  shown  by  both  peaks  are  well  above  the 
factor  of  two  predicted  and  found  in  prior  studies  involving  extended  scattering  media.  Such  greater- 
than-expected  backscattering  enhancement  remains  unexplained  at  the  current  time. 


<j»  (degrees) 


/,  ,2  \  / /i  |2  \ 

Figure  2.  Mean-square  scattered-field  ratio  ( k t/s(R,  (p)  y/ (/?,  (p)  )  vs.  scattering  angle  <pfor 

\  1 1  \  /  ave 

four  different  average  scatterer  spacings:  k0s  =  11. 7  (blue),  5.5  (red),  3.2  (green),  and  2.5  (purple). 
Here  scatterers  were  placed  in  a  cubical  aggregation  with  5  wavelength  edges.  The  ensemble 
average  was  computed  from  256  to  1024  Monte-Carlo  trials.  The  acoustic  coherent  backscatter 
enhancement  shown  here  for  kos  =  3.2  and  2.5  is  above  that  expected  from  prior  optical  and  acoustic 

studies  involving  extended  scattering  media. 
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When  the  comers  and  edges  of  the  cube  are  progressively  removed,  the  ACBE  peak  amplitude 
decreases  and  the  peak  width  mildly  increases  as  the  aggregation  geometry  becomes  progressively 

more  spherical.  Figure  3  shows  such  results  for  kQs  =  2.5  and  k0Jof=  3.5.  The  vertical  axis  is 

/,  ,2  \  If  |  ,2\ 

\f/s{R,(p)  y/s(R,<p)  and  the  horizontal  axis  is  the  far  field  scattering  angle  <p.  The  peak 

\  /  /  \  /  ave 

amplitude  decreases  as  the  aggregation  becomes  more  spherical.  The  amplitude  for  the  spherical 
aggregation  (grey  curve)  is  closest  to  the  enhancement  expected  from  prior  studies. 


\ D  (degrees) 
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Figure  3.  Mean  square  scattered-field  ratio  (  y/fR,  (p)  y/ (/?,  cp)  vs.  scattering  angle  (pfor 

\  1 1  \  i  ave 

an  average  scatter er  spacing  k„s  =  2.5.  Here  scatter ers  are  initially  placed  in  cubical  aggregation 
with  a  5-wavelength  edge  (blue),  and  the  corners  and  edges  of  the  scattering  aggregation  are 
progressively  removed  by  dropping  scatters  beyond  progressively  smaller  distances  from  the  center 
of  the  aggregation.  Results  are  shown  when  scatterers  are  dropped  at  distances  greater  than  42, 
(red)  3. 72,  (green)  3.32  (purple)  32  (light  blue),  2. 72  (orange)  and  2.52  (gray).  The  amplitude  of  the 
peak  at  (p  =  0  decreases  as  the  corners  and  edges  of  the  scattering  aggregation  are  removed. 

As  a  further  attempt  to  determine  if  these  predictions  are  correct,  the  simulation  geometry  was  adjusted 
to  match  a  combination  of  parameters  from  two  CBE  experiments  and  then  compared  to  relevant 
theoretical  results.  The  scattering  media  parameters  are  obtained  from  an  optics  experiment3  while  the 
simulation  geometry  is  that  of  an  ultrasound  experiment.4  Both  of  these  experiments  need  to  be 
considered  in  order  to  accommodate  the  limitations  associated  with  the  current  simulations.  The  optics 
experiment3  utilized  a  small  wavelength  and  a  large  number  (trillions)  of  nominally-spherical  scatterers 
in  a  three  dimensional  geometry.  However,  the  current  simulation  cannot  be  completed  for  such  a  large 
number  of  scatterers.  The  ultrasound  experiment4  utilized  a  few  thousand  scatterers  -  a  number  that 
can  be  managed  by  the  current  simulations  -  but  they  were  two-dimensional  scatterers  (vertical 
metallic  rods).  Wave  scattering  from  such  long  cylinders  is  different  than  that  from  the  omni¬ 
directional  point  scatterers  in  the  current  simulations.  So,  in  order  to  bridge  the  gap  between  optics  and 
acoustics,  the  relevant  scattering-medium  length-scale  ratios  of  the  comparison  simulations  are 
obtained  from  the  optics  experiment  while  the  scattering-aggregation  and  data-collection  geometries 
are  that  of  the  ultrasound  experiments. 
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These  simulations  were  conducted  using  the  geometry  shown  in  Fig.  1  with  the  following  parameters: 
=  0.232,  kos  =  3.13,  ko R  =  131,  k® L  =  235,  k<Ji  =  293,  and  k$w  =  88.  The  solid  line  is  the 

simulation  result  and  the  dashed  line  is  the  theoretical  CBE  result  based  on  the  mean-free  path  in  the 
scattering  medium.  The  results  shown  in  this  figure  required:  (/)  fitting  the  mean-field  intensity  decay 
inside  the  scattering  aggregation  of  the  simulations  to  determine  the  mean  free  path,  (ii)  matching  the 
scattering  efficiency  and  spacing  of  the  polystyrene  spheres  in  the  optics  experiment  to  the  simulation's 
point  scatterers,  and  (iii)  accounting  for  the  finite  angular  resolution  of  the  linear  receiving  array. 
Overall  the  match  shown  in  Fig.  3  is  good.  The  minor  differences  between  the  two  curves  are  likely 
related  to  uncertainty  in  the  fitting  for  the  acoustic  mean  free  path  over  the  finite  extent  ( w  in  Fig.  1)  of 
the  simulated  scattering  aggregation. 


(j)  (degrees) 


Figure  4.  Comparison  of  simulated  and  theoretical  results  for  ACRE  using  a  composite  selection  of 
experimental  parameters.  The  vertical  axis  is  B($)/[B]ave  and  the  horizontal  axis  is  the  beam  steering 
angle  (j>  shown  in  Fig.  1.  The  solid  line  is  the  simulation  result.  The  dashed  line  is  an  evaluation  of 
the  relevant  theory2.  The  minor  differences  between  the  two  curves  are  likely  related  to  uncertainty 
in  determining  the  acoustic  mean  free  path,  a  parameter  in  the  theory  that  must  be  obtained  by 
fitting  the  mean-field  intensity  decay  inside  the  scattering  aggregation  of  the  simulations. 


Overall,  these  checks  of  the  simulations  have  not  yet  found  any  deficiencies  or  errors.  Thus,  the 
greater-than-expected  ACBE  peak  values  may  be  genuine.  However,  an  explanation  for  these  high 
peak  values  has  not  yet  been  found,  but  the  search  continues  and  currently  centers  on  possible  coherent 
scattering  and  aggregation  shape  effects  not  included  in  CBE  theory.  A  journal  manuscript  will  be 
prepared  when  this  explanation  has  been  found. 

IMPACT/APPLICATION 

In  broad  terms,  this  project  ultimately  seeks  to  determine  if  and  how  ACBE  might  be  exploited  for 
active  sonar  applications.  In  particular,  if  successful,  it  should  prove  useful  for  remote  classification, 
because  a  large  sonar  return  from  a  single  large  scatterer  will  likely  not  display  ACBE  while  a 
similarly  large  sonar  return  from  an  aggregation  of  many  small  scatterers  may  display  ACBE.  Thus, 
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this  research  effort  may  eventually  impact  how  active  sonar  signals  are  processed  and  displayed  for 
tactical  decision-making  related  to  classification. 

TRANSITIONS 

The  results  of  this  research  effort  should  aid  in  the  design  of  active  sonar  signal  processors  for  tactical 
decision  aids.  However,  at  this  time  no  direct  transition  links  have  been  established  with  more  applied 
research  programs.  Once  the  current  simulation  capability  is  more  firmly  established  and  validated, 
and  promising  results  have  been  obtained,  a  transition  path  through  NRL  or  one  of  the  Navy's  Warfare 
Centers  will  be  sought. 

RELATED  PROJECTS 

This  project  is  related  the  other  projects  funded  under  ONR’s  2010  basic  research  challenge  program. 
In  particular,  the  work  by  Prof.  Feuillade  in  Chile  is  most  closely  related. 
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